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ABSTRACT. The results of a search for distant clusters of galaxies performed 
using the I-band data obtained by the ESO Imaging Survey (EIS) are pre- 
sented. Cluster candidates are identified using a matched filter algorithm, 
that provides not only an objective detection criterion, but also the means 
to estimate the cluster redshift and richness. A preliminary sample of distant 
clusters has been obtained, containing 252 cluster candidates with estimated 
redshift in the interval 0.2 < z < 1.3 (median redshift z mc( j ~ 0.4) over an 
area of approximately 14 square degrees. The adopted selection criteria for 
the inclusion of cluster candidates in this sample has been in general con- 
servative, as the primary concern has been the reliability of the candidates 
rather than the completeness of the sample. 



1 Introduction 

One of the primary goals of the re- 
cently completed ESO Imaging Sur- 
vey (EIS; Renzini & da Costa 1997) 
has been the preparation of a sample 
of optically-selected clusters of galax- 
ies over an extended redshift base- 
line. Such a sample could be used for 
many different studies, ranging from 
the evolution of the galaxy popula- 
tion in clusters, to the evolution of 
the abundance of galaxy clusters, a 
powerful discriminant of cosmolog- 
ical models. Until recently, only a 
handful of clusters were known at 
redshifts z > 0.5; visual searches for 
high redshift clusters were conducted 
by Gunn etal. (1986) and Couch 
etal. (1991), but their samples are 
severely incomplete beyond z ~ 0.5; 
at higher redshifts targeted obser- 
vations in fields containing known 
radio-galaxies and QSOs have pro- 
duced a handful of cluster identifi- 
cations (e.g. Dickinson 1995; Fran- 
cis etal. 1996; Pascarelle etal. 1996; 
Deltorn etal. 1997). The first objec- 
tive search for distant clusters was 
conducted by Postman etal. (1996; 
hereafter P96) using the 5-m tele- 
scope of the Palomar Observatory. In 
that survey 10 of the 79 cluster can- 
didates that were identified have es- 
timated redshift > 0.7. Further ev- 
idence for the existence of clusters 
at high redshift has been obtained 
from X-ray (e.g. Gioia & Luppino 
1994; Henry et at 1997; Rosati et al. 



1998), optical (e.g. Connolly etal. 
1996; Zaritsky etal. 1997) and in- 
frared (Stanford etal. 1997) searches. 
However, the existing samples are 
still small, and their selection effects 
largely unknown. 

In this paper we briefly describe 
a large sample of optically-selected 
cluster candidates, obtained apply- 
ing a matched filter algorithm to the 
catalogs of objects detected on sin- 
gle 150 sec. EIS V- and I-band expo- 
sures. The sample includes 252 can- 
didates, identified over an area of 
14.4 deg 2 . Further details on this 
work can be found in Olsen et al. 
(1998a,b; hereafter 098a,b), and in 
Scodeggio etal. (1998). 

2 Observations and Data Re- 
duction 

EIS observations were conducted 
using the EMMI camera on the 
ESO 3.5m New Technology Tele- 
scope. Observations were carried out 
over four pre-selected patches of the 
sky, spanning a wide range in right 
ascension, and cover a total area of 
17 deg 2 in I-band. Only small sub- 
areas were covered also by V- and B- 
band observations (3.5 and 1.5 deg 2 , 
respectively) . For this reason we con- 
centrate here our attention mostly 
on the I-band data. EIS observations 
consist of a sequence of 150 sec expo- 
sures; however each point of a patch 
is imaged twice, for a total integra- 
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tion time of 300 sec, using two frames 
shifted by half an EMMI-frame both 
in right ascension and declination. 
The easiest way of visualizing the 
global geometry of this mosaic of 
frames is to consider two indepen- 
dent sets of them, forming contigu- 
ous grids (in the following referred to 
as odd and even frames), superposed 
and shifted by half a frame both in 
right ascension and declination. 

The data reduction is carried 
out automatically through the EIS 
pipeline, described in Nonino et al. 
(1998; hereafter N98). This pipeline 
produces both coadded images and 
fully corrected single frames, us- 
ing the astrometric and photomet- 
ric solution derived from the global 
data reduction process. The astro- 
metric solution is found relative to 
the USNO-A1 catalog. The inter- 
nal accuracy of the astrometric so- 
lution is better than 0.03 arcsec, al- 
though the absolute calibration suf- 
fers from the random and systematic 
errors of the reference catalog. The 
photometric calibration is done in a 
two step procedure, first bringing all 
the frames to a common photomet- 
ric zero-point, taking advantage of 
the overlap between the frames, then 
making an absolute calibration based 
on external data. The internal accu- 
racy of the photometric calibration 
is jCO.005 mag. The current abso- 
lute calibration uncertainty is <^0.2 
mag. Further details can be found in 
N98. 

Among the various object cata- 
logs produced by the EIS pipeline, 
we use in this work the so-called odd 
and even catalogs, which are single 
entry catalogs listing all objects de- 
tected in the odd or even frames. 
The reliability and completeness of 
these catalogs are discussed in N98. 
Based on that analysis, it was esti- 
mated that the single-frame odd and 
even I-band catalogs are ~95% com- 
plete to i" = 23.0; with a differen- 
tial completeness at this magnitude 
of 80% (the V-band ones have sim- 
ilar completeness to V = 24.0). At 
this limiting magnitude the contam- 
ination from spurious objects is es- 
timated to be approximately 20%, 
with total contamination of the cata- 
log of ~5%. The object classification 
was shown to be reliable to / Ri 21. 
Below that limit the object classifica- 
tion is not reliable any more, and all 



detected objects fainter than I = 21 
are taken to be galaxies. Already at 
this magnitude the fraction of stars is 
found to be ~25% of the total num- 
ber of objects, and taking into ac- 
count the steep rise of the galaxy 
number counts faintward than / = 
21, the contamination of the galaxy 
catalog by stars can be considered 
negligible. 

3 Cluster detection 

Candidate galaxy clusters are ob- 
jectively identified using a procedure 
based on the matched filter algo- 
rithm described by P96. The details 
of this EIS cluster finding pipeline 
are discussed by 098a, and only a 
brief summary is given here. The fil- 
ter is derived from an approximate 
maximum likelihood estimator, ob- 
tained from a model of the spatial 
and luminosity distribution of galax- 
ies within a cluster. This distribu- 
tion, as a function of magnitude m 
and radial distance from the cluster 
center r c , is represented as 

C(r, m) = A cl P(r/r c )(f>(m - m*) (1) 

where P is the cluster projected ra- 
dial profile, 4> is the cluster lumi- 
nosity function, and A c i measures 
the cluster richness. The parameters 
m* and r c are the apparent mag- 
nitude corresponding to the char- 
acteristic luminosity of the cluster 
galaxies, and the projected value 
of the cluster characteristic scale 
length. From this model one can rep- 
resent the observed galaxy distribu- 
tion DasD(r,ra) = b(m) + C(r, m) 
where b(m) is the background num- 
ber counts, and write an approxi- 
mate likelihood C of having a cluster 
at a given position as 

ln£~ / C(r ' m) D(r, to) d 2 r dm(2) 
J b ( m ) 

The functional dependence of C 
on the redshift-dependent parame- 
ters m* and r c provides the mean 
of estimating the cluster candidate 
rcdshift. By evaluating the func- 
tion C for each element of a two- 
dimensional array a filtered image 
of the input galaxy catalog is cre- 
ated (the "Likelihood map"). This 
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procedure is repeated a number of 
times, with the parameters m* and 
r c tuned for different cluster red- 
shift values. Significant peaks in the 
likelihood maps are identified inde- 
pendently in each map, using SEx- 
tractor. A minimum size criterion is 
adopted, whereby only peaks where 
more than 7V m i n pixels have values 
above the detection threshold cr^ ct 
are considered as potential detec- 
tions. The value of iV m ; n is set to cor- 
respond to the area of a circle with 
radius lr c , while the value of tTjet is 
kept constant at 2. These parameters 
were optimized using simulations de- 
scribed in 098a. 

The lists of peaks identified in 
the various Likelihood maps are then 
compared, and peaks detected at 
more than one filter redshift are asso- 
ciated on the basis of positional coin- 
cidence. From this association, like- 
lihood versus z curves are created, 
and the redshift and richness esti- 
mates for each candidate are derived 
locating the peak of the correspond- 
ing likelihood versus z curve. Two 
richness parameters are derived for 
each candidate, following P96. The 
first is obtained from the matched 
filter procedure itself, using the pa- 
rameter A c i, while a second indepen- 
dent richness estimate is obtained 
to reproduce more closely the con- 
ventional Abell richness parameter. 
The final cluster candidates sample is 
then composed of those objects that 
persist for at least three filter red- 
shift tunings, and that have a rich- 
ness estimate A c i > 30. Simulations 
have shown these to be quite con- 
servative selection criteria, as noise 
peaks are most likely to be associ- 
ated with small values of the inferred 
richness, and to be detected only for 
one (or few) particular redshift tun- 
ing of the matched filter (see 098a 
for details). 

4 Results 

The total sample of EIS cluster 
candidates, detected in 1-band over 
the four EIS patches, consists of 252 
objects with estimated redshift in the 
range 0.2 < z < 1.3. The total 
area over which the search was car- 
ried out is of 14.4 deg 2 , and there- 
fore the inferred density of candi- 
dates is of 17.5 deg -2 . This value 



is slightly larger than that found 
by P96, and could result from the 
somewhat deeper magnitude limit 
achieved by EIS, with respect to the 
P96 data. The distribution of esti- 
mated rcdshifts for this total sample 
is shown in Figure [l]. The median es- 
timated redshift is z ~ 0.4, and coin- 
cides with the value for the P96 sam- 
ple. Note however that the two red- 
shift distributions differ somewhat, 
with the EIS sample showing a more 
extended tail beyond z ~ 0.6 (upper 
panel of Figure 0) . A similar compar- 
ison can be carried out also with the 
sample of X-ray selected clusters of 
Rosati etal. (1998). Using the cluster 
confirmations and redshift measure- 
ments obtained up to December 1997 
for that sample (Rosati, private com- 
munication) , we can see once again a 
general agreement (the median red- 
shift of the clusters in the Rosati 
etal. sample is also z ~ 0.4), and 
a more extended tail of high red- 
shift objects in the EIS sample (lower 
panel of Figure ^). Spectroscopic ob- 
servations are now necessary to de- 
termine what fraction of EIS cluster 
candidates are real objects, and to 
determine more accurately their red- 
shift. 

For applications of cosmological 
interest, it is also interesting to con- 
sider the distribution of rich (and 
presumably most massive) clusters. 
This is presented in Figure |l] as 
the shaded portion of the histogram. 
From this distribution it is possible 
to deduce that the comoving density 
of cluster candidates is changing only 
very mildly as a function of redshift. 

Over an area of 2 square deg. 
the search for cluster candidates has 
been performed also using V-band 
data, independently from the search 
in I-band. Afterwards the V and I- 
band candidate lists were compared, 
and V-band detections associated to 
I-band ones. Of the 19 candidates 
with estimated redshift z < 0.5, 17 
(~ 90%) are also detected in V-band; 
for those with estimated redshift z > 
0.5, only 4 out of 16 are detected in 
V-band. This result is not surprising 
since the V-band data are relatively 
shallower the I-band ones. Using the 
rule of thumb that the data should 
reach at least one magnitude fainter 
than m* at the redshift of a given 
cluster to allow for its detection, one 
can translate the galaxy catalog lim- 
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Figure 1. The estimated redshift distribution for the total sample of 252 
EIS cluster candidates, and for the subset of 111 candidates with richness 
estimate > 40 (shaded portion of the histogram). 
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Figure 2. Comparison of the estimated redshift distributions for the total EIS 
sample and for the sample of Postman etal. (upper panel); and comparison 
between the estimated redshift distribution of the EIS sample and the mea- 
sured redshift distribution of the Rosati etal. sample (lower panel). The EIS 
sample is scaled down to the same area covered by the Postman et al. survey, 
and to the same total number of clusters as in the Rosati et al. survey. 



iting magnitudes adopted in the clus- 
ter search into limiting redshifts for 
cluster detection. The I-band limit 
of / = 23.0 then translates into a 
limiting redshift between ~ 1.0 (no- 
evolution model) and ~ 1.3 (pas- 
sive evolution model), whereas the 
V-band limit of V = 24.0 translates 
into a limiting redshift between ~ 
0.55 (no-evolution) and ~ 0.7 (pas- 
sive evolution). 

The availability of data in two 
passbands provides also an alterna- 
tive way of confirming cluster candi- 
dates and their estimated redshifts, 



based on the detection of the se- 
quence of cluster early-type galax- 
ies in a Color-Magnitude (CM) di- 
agram. In order to investigate this 
possibility, a V-I vs I CM diagram 
was produced for each cluster can- 
didate, showing all galaxies within 
a radius of 0.75 h _1 Mpc (Ho = 
75 km s _1 /Mpc) from the nominal 
cluster center. Fig. ^ shows four ex- 
amples of such diagrams, illustrating 
cases with estimated redshift in the 
range 0.2 < z < 0.6. At low redshift, 
the sequence of early-type galaxies is 
clearly visible, but at z ^>0.5 the ev- 
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Figure 3. Color-magnitude diagrams observed for four cluster candidates 
with estimated redshifts in range of z = 0.2 to z = 0.6, as indicated in each 
panel. Also shown are the values of m| and the expected colors of typical 
ellipticals at these redshifts. 



idence for a CM relation is, in most 
cases, less compelling. In total, out of 
35 cluster candidates in the region of 
overlap of the V- and I-band images, 
there are 19 with evidence for a CM 
relation, with estimated redshifts ex- 
tending out to z ^C0.6. Furthermore, 
the redshift estimates based on the 
color and on the matched filter seem 
to agree, in most cases, within 0.1. 
However, there are at least four cases 
where there is a strong suggestion 
that the matched filter has overesti- 
mated the redshift. 
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